Mesoscopic p-wave superconductor near the phase transition temperature 
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We study the finite-size and boundary effects on a p-wave superconductor in a mesoscopic rectan- 
gular sample using Ginzburg-Landau (GL) and quasi-classical (QC) Green's function theory. Except 
for a square sample with parameters far away from the isotropic weak-coupling limit, the ground 
state near the critical temperature always prefers a time-reversal symmetric state, where the order 
parameter can be represented by a real vector. For large aspect ratio, this vector is parallel to 
the long side of the rectangle. Within a critical aspect ratio, it has instead a vortex-like structure, 
vanishing at the sample center. 

PACS numbers: 74.78.Na, 74.20.De, 74.20.Rp 



Studies of multicomponent superfluids and supercon- 
ductors have excited many for decades because of the 
diversity of textures, complex vortex structures and col- 
lective modes. The superfiuid 3 He with spin-triplet or- 
der parameter [HQ is a well-established example. Many 
studies also show that superconductors with multicom- 
ponent order parameters can also be found in, for exam- 
ple, UPt 3 [111] and Sr 2 Ru0 4 [1,11]. Recently, studies of 
multicomponent superconductor in a confined geometry 
draw much attention due to advancements in nanofabri- 
cation. Experiments claimed to find half-quantum vor- 
tices and the Little-Parks effect Q in Sr 2 RuC>4 quan- 
tum ring. Surfaces are expected to have non-trivial ef- 
fects on such superconductors. Some theoretical works 
show that surface currents are present in broken time- 
reversal symmetric superconductors. [9MTTj] In consider- 
ing Ru inclusions, Sigrist and his collaborators [12] have 
shown that a time-reversal symmetry state can be favored 
near the interface between Ru and Sr 2 Ru04 due to the 
boundary conditions. In our previous work [l3| . consid- 
ering a thin circular disk with smooth boundaries and 
applying Ginzburg-Landau (GL) theory, we have shown 
that a two-component p-wave superconductor can exhibit 
multiple phase transitions in a confined geometry. At 
zero magnetic field, the superconducting transition from 
the normal state was found to be always first to a time- 
reversal symmetric state (with one exception which we 
discussed in more detail below), even though the bulk 
free energy may favor a broken time-reversal symmetry 
state, which can exist at a lower temperature. This time- 
reversal symmetric state has a vortex-like structure, with 
order parameter vanishing at the center of the disk. We 
have also argued there that these features are general, do 
not rely on the GL approximation and should exist for 
also general geometries. [l4| 

In this paper, we shall show that this is indeed true 
by considering rectangular and square samples. For rect- 
angular samples with large aspect ratios, we show that 
the phase transitions are first to a state with order pa- 
rameter being a real vector parallel to the long side of 
the sample. For smaller aspect ratios, the state near the 
transition temperature is again a time-reversal symmet- 
ric state with a vortex at the center, much like what we 



found for the circular disk (again with an exception to 
be mentioned below). We show these results using both 
GL and quasiclassical (QC) theories. The results from 
these two approaches are qualitatively similar except for 
the critical sizes and aspect ratios obtained. 

We shall thus consider a superconductor where its or- 
bital part is given by ff = r\ x x + rj y y. We shall consider 
the dependences of r\ x and r\ y on the coordinates x, y, 
assuming that they are constant along the z direction. 
We assume the length of the sample in x direction is L 
and the width in y direction is W, and these surfaces 
are smooth. The effects of rough boundary have been 
discussed in Ref.[l3[] for the circular disk. We shall also 
limit ourselves to zero external magnetic fields. Near the 
transition temperature, the magnetic field generated by 
the supercurrent is also negligible, hence the vector po- 
tential can always be ignored. 

First, we study this system via GL theory. The GL 
free energy density per unit area for the bulk, J-),, can be 
written as 



T h = a(rf* ■ ff) + 



(1) 



where a = a'(t - 1) with a' > 0, t = T/T c ° is the ra- 
tio of the temperature T relative to the bulk transition 
temperature T®, and . . . represents terms higher power 
in the order parameter which are irrelevant below since 
we are interested only in the physics at the (modified) 
transition temperature T c . In the presence of spatial 
variations, there is an additional contribution to the free 
energy given by 



T g = K^djrnWjmy + K 2 (d jVj )(d im )* 

+K 3 (d j r ]l ){d l r 1j )* + K 4 [(d xVx ) 2 + (dyriy) 2 } 



(2) 



where repeated indices j, I in the first three terms are 
summed over x, y, and the last term describes crys- 
tal anisotropy. [1 5j Within weak coupling approximation, 
particle-hole symmetry, and for an isotropic Fermi sur- 
face, K\ = K2 = K3 > and K4 — 0, but wc shall treat 
these coefficients as general parameters. 

The GL equations need to be accompanied by bound- 
ary conditions. The perpendicular component of the or- 
der parameter at the surface should vanish [17) . Thus, 



2 



for a point at the surface where the normal is h, h-fj = 0. 
For a smooth surface, the parallel component m should 
have vanishing normal gradient [13] • That is, at the sur- 
face, {h ■ V)f7|| = 0. 

GL equations for r\ x ^ y can be obtained by the variation 
principle. Near the critical temperature, we can linearize 
these equations. The easiest way to match the bound- 
ary conditions is to superimpose the Fourier components. 
Written in matrix form, the GL equations for the Fourier 
component q become 



K iq 2 + K 234 q x 
K 23 q x q y 



K 23 q x q y 

K ig 2 + K 23i q 2 y 



Vx.q 



= a'(l -t) 



Vx,q 



(3) 



Here, q is the wavevector, q 2 = q 2 + q 2 , K 23 = K 2 + K 3 , 
and K 234 — K 2 + K 3 + if 4. 

It is easy to find that, if q x — or q y — 0, Eq.Q 
decouples. We obtain either (i) r\ x ^ 7^ with r\ y ^ = 
or (ii) r\ y ^ 7^ with i] x ^ — 0. We call these solutions 
as A phases. In case (i), we have two possibilities. One 
is q = x and the critical temperature is determined by 
a'(l — t) = K\ 234 q 2 • Because r) x (x) is independent of y, 
the possible solutions are rj x = X sin , satisfying the 
boundary conditions, rj x = 0, at x = and L. Here X is 
a constant and m is an integer. The best choice is to = 1 
and the critical temperature is a'(l — t) = K\ 234 {jt j V) 2 ■ 
We call this the A\ phase. The other is q — y. The 
order parameter rj x (y) is independent of x. Thus it is 
not possible to satisfy the boundary conditions at x = 
and L. In case (ii), the best solution is rj y oc sin which 
is just the solution in case (i) with x ■(-> y. The critical 
temperature is determined by a'(l —t) = Ki 23 4(tt /W) 2 . 
We call this the A 2 phase. 

If both q x and q v are non-zero, both r\ x ^, r\ y ^ are finite. 
We define this kind of solution as B phase. To simplify 
the calculations, we ignore crystal anisotropy for the mo- 
ment, and set K4 = 0. From Eq.((31), we find the smallest 
eigenvalue is K\q 2 (for K 23 > 0). To have the normal 
component to the surfaces at x = and L to vanish, r\ x 
must have the factor sin(m7rx/L), where m is an integer. 
Because the boundary conditions drj x /dy — at y = 
and W, rj x should be proportional to cos(mry/W), with 
n also an integer. We can use the same arguments for 
7] y . Therefore 



rj x = X sin cos 
rjy = Y cos THIS. s j n 



mry 

~W' 



(4) 



The critical temperature is highest for m = n = 1, and 
thus determined by a'(l - t) = K^tt/L) 2 + (n/W) 2 }. 
In order to satisfy Eq.©, we need (n/L)X + (n/W)Y = 
0, which means X and Y are relatively real and have 
opposite signs. 

Comparing the transition temperatures of the A and B 
phases, we find that the system prefers the A\ phase for 
L ^> W. For L ~ W, it prefers the B phase. We define 
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FIG. 1: GL phase diagram at the transition temperature. 
K4. = 0. (a) K 23 /K! > 1 (b) K 2 3/K! < 1. 



the aspect ratio of the sample as p = L/W. The critical 
aspect ratio separating these two phases is 



Pc = ^/Kzz/Kx. 



(5) 



With the same reasoning, the system is in A 2 phase for 
W ^> L but prefers the B phase if p is larger than 
p^ 1 = y/ K\ I K 23 . The phase diagram is shown in 
FigUJa). The cartoon pictures under the ruler are used 
to specify main characteristic of the order parameter in 
the corresponding phases. The solution for the middle 
case (B phase) is qualitatively same as the circular disk 
in [13j. Because the relatively real coefficients in Eq.@, 
we can represent rj by a real vector, as done in the inset 
of FigJ^fc). It is clear that the order parameter forms a 
vortex-like structure, and vanishes at the sample center. 

When crystal anisotropy is included, the critical ratio 
becomes 



'-^1234^234 - ^4(^23 + -?^234) 



(6) 



y K\K\ 234 

It reduces to Eq.© for K4 = 0. For small K 4 , the critical 
ratio is \/2 — K 4 / K\ 23 within the weak-coupling limit. It 
shows that the crystal anisotropy for K4 > (<)0 stabi- 
lizes(destabilizes) the order parameter with the direction 
parallel to the long side of the sample. The phase dia- 
gram is similar to FigQTa) except smaller(larger) region 
for vortex state. We ignore crystal isotropy in the follow- 
ing. 

With decreasing K 23 /K 1: the stability region for B 
phase becomes narrower. This phase diagram (FigQTa)) 
will change qualitatively if K\ > ^"23 The B phase is 
never stable (at T c ), and the system is in the Ai phase 
if p > 1, and in the A 2 phase if p < 1. The square 
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sample p — 1 forms a special case, where the system 
still has C4 symmetry in real space, and the Ai and A 2 
phases are therefore degenerate. One can combine these 
two solutions with a phase difference. As a result, if the 
higher order terms in eq ([T]) for the bulk free energy prefer 
a timc-reversal-symmetry-broken state, then the system 
would enter such a state directly at T c . Therefore, the 
phase diagram becomes FigQJb) f° r -^23 < Ki, where 
the ground state at p = 1 should break the time-reversal 
symmetry. [l8j 

Our results for the square here provide further under- 
standing of those we obtained earlier for the circular disk 
There, we found that, for -ftTi.2,3 near the weak- 
coupling values (Ki = K% = K 3 ) the phase transition 
from the normal state is always to the state (named n = 1 
there) which preserves time-reversal but with a vortex at 
the center. That phase obviously has the same qualita- 
tive behavior as our B phase here. For sufficiently small 
K23/K1, we found that the system can enter a broken 
time-reversal symmetry state directly. We find the same 
results here though the critical value for K23/K1 obvi- 
ously can depend on the geometry. 

In order to check the validity of the phase diagram from 
GL theory, we employ quasi-classical (QC) Green's func- 
tion theory. For simplification, we focus on the isotropic 
and weak-coupling case. As shown in Ref.[l3|, we have, 
after linearizing in the order parameter, 



(2ie n + ivfp ■ V)/ = 2i7r(sgne„)A, 



(7) 



where f(p,e n ,f) and A(p, r) describe separately the off- 
diagonal parts of the QC propagator and pairing func- 
tion, p is the momentum direction, e„ is Matsubara fre- 
quency, and Vf is the Fermi velocity. With pairing inter- 
action written as V\p • p', the gap equation reads 

A(p, r) = N(0)TVi J2<(P- P')f(P'' e») > (8) 



where the angular bracket denotes angular average over 
p' and N(0) is the density of states at the Fermi level. 
For our square geometry and assuming smooth surfaces, 
we have the boundary conditions f{9) = f(n— 0) at x=0 
and L and f(6) = f(-9) at y=0 and W. Here 9 is the 
angle between p and x. 

Before solving the case in a confined rectangle, we like 
to mention the connection between GL theory and QC 
theory for the bulk. To zeroth order in gradient, one finds 



2e„ 



(9) 



which defines the bulk transition temperature T®. The 
first order for / is odd in e„ and will not contribute to 
the gap equation. In the second order, we recover the GL 
theory with 



c ^41 



< cos 



' sm 



>, 



(10) 



and similar expressions for if 2,3, with K\ = Ko = K 3 . 
Our equations are consistent with those in Ref. [H H, Hq 
For the A\ phase, we shall show that we can have a 
self-consistent solution in QC theory with 



A(p, r) = X sin — cos 9, 



(11) 



the order parameter suggested by the GL theory. One 
finds that / is independent of y. With the ansatz 

f(0,e n ;x) =Ci(0,e„)sin— + C 2 (0, e„) cos — (12) 

satisfying the boundary conditions, solving for Ci^ifi, e„) 
via eq.JTJ) and using ©, we find 



, T° ^ (^f) 2 cos 4 

ln-f = 2^T c J2 



' " ^1 I r-t , ( ^7 — COS 0) 2 

e „=-oo 4|e»| 3 [l+ \ £j 



(13) 



For large L, one can replace the bracket in the denom- 
inator by 1, the LHS by (1 — t), recovering the GL re- 
sult using Eq. fTU)) . Hence we see that, beyond GL, one 
needs simply to include extra factors in the denominator 
of Eq. ([T3"]) and include the In on the LHS. 

Now, we consider the B phase. The order parameter is 

. , „ . irx Try irx , iry . 

A(p, r) = X sm — cos — cos 9 + Y cos — sm — sm 9. 

lj W Li W 

(14) 

From previous experience, it suggests that the solution 
has the following form 

/ = e„) sin S cos 3 + C 2 (9, e n ) cos ^ sin 

C 3 (9, e n ) cos % cos §f + C 4 (0, e n ) sin ^ sin at. 

To simplify writing, let A = irvf/L, B = nVf/W. Again 
solving for / from (J7|, we obtain the following coupled 
linear equations in X,Y: 



X\Jl - 2^T C ^( 



[(ci +c 2 )cos 2 9}X + c 3 Y 



Yln- 



\e n \D 

2 , 



2irT c 



,[( Cl + c 2 )sin" 9}Y + c 3 X, 
\e n \D 



(16) 



(17) 



Here c a = (A 2 cos 2 9 + B 2 sin 2 6/)/(4e 2 l ), c 2 = (A 2 cos 2 (9- 
B 2 sin 2 6») 2 /(4e 2 ) 2 , c 3 = (AB sin 2 9 cos 2 9)/(2e 2 l ), and 
D = 1 + 2c\ + c 2 . The critical temperature of the B 
phase is determined by the point which allows non-trivial 
X and Y. We note that if one keeps only the lowest or- 
ders in A 2 and B 2 , D — > 1, and replaces the ln's on LHS 
by (1 — t), then Eq. (fP6]) and Eq. (fT7|) recover the corre- 
sponding equations ([3]) in GL theory. 

In Figf2][a), we compare the critical temperatures for 
different sizes of square samples between GL and Q C 
theories. We use the coherence length £ = \jKx13j01l as 
the unit for length (£ = 0.199909w//T c ° for QC). In GL 
theory, (1 — t), the relative suppression of critical tem- 
perature, is inversely proportional to the square of the 
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FIG. 2: (a) Critical temperatures from GL and QC theories 
for square samples, (b) Critical temperatures in QC theory 
for different p's. As the size of the system is reduced, the 
ground state near the critical temperature can change from A 
to B. The possible first order phase transition is indicated in 
the region between dash lines, (c) Phase diagram for different 
sizes. Lines with squares or circles correspond to the second 
order phase transition; a first order transition line lies within 
each dashed region. Inset: Order parameter for a square sam- 
ple. Note that all results are for K\ = K2 = K3. £ is the 
coherence length. 



length scale of the system. Therefore it is more conve- 
nient to set the vertical axes of the phase diagram to be 
(£/!/) 2 . We obtain the straight line with crosses for the 
critical temperature of A phase and the line with pluses 
for that of B phase. It shows that the system prefers the 
B phase. On the other hand, we also present the crit- 
ical temperatures calculated from QC theory. The line 
with squares is for the B phase and the line with circles 
is for the A phase. As expected, it shows that the re- 
sults from QC theory are consistent with those from GL 
theory near t = 1 which is corresponding to the bulk sys- 
tem —> 0). In addition to the fact that the B phase 
is still preferred, we see that the critical temperature is 
more suppressed than GL as the size of the system is 



smaller. [19( 

As the aspect ratio p becomes larger than 1, GL the- 
ory shows that there is a phase transition from the B 
phase to t he A ph ase as p increases beyond the criti- 
cal ratio y / K 23 /K 1 (Fig.l(a)). In Fig.2(b), we find that 
this critical ratio (\/2 here) still applies for large systems. 
However, we find that the B phase occupies a slightly 
larger p region when the size decreases. An example is 
in Fig. 2(b), where we show that the system processes a 
phase transition from the A phase (thick black line) to 
the B phase for p = 1.6 (thin red line) around t = 0.67 
as the size of the system becomes smaller. As the aspect 
ratio decreases further, the situation becomes more com- 
plicated because the phase boundary for A phase is not 
a monotonic function in t. From the shape of the curve, 
we expect that we have a first order phase transition at 
lower temperatures (cf. 20] ) . We shall leave the detailed 
investigation of this first order phase transition for the 
future. For illustrative purposes, we indicate the result- 
ing phase diagram by postulating a first order transition 
line somewhere between the two dashed lines in Fig. 2(b). 
Our obtained phase diagram is summarized in Fig. 2(c). 
Here, instead of using (£,/L) 2 as the vertical axes, we 
use the area of the sample. Fig. 2(c) shows that for very 
small sample within the QC limit, the system is normal 
and not superconducting. As the size of the system is 
larger, one can get the A phase via a first order phase 
transition for sufficiently large or small aspect ratios, or 
otherwise the B phase via a second order phase transi- 
tion. The B phase with the aspect ratio larger than \/2 
or smaller than 1 / \[2 can have another phase transition 
into A phase. For very large samples, we recover the crit- 
ical ratios V2 and l/v 7 ^, consistent with the results from 
GL theory. For the square samples, the stable phase is 
the B phase, with a vortex-like structure at the center, 
similar to the n = 1 state in [l3j obtained in the disk 
geometry. 

In conclusion, we studied a two-component p-wave su- 
perconductor in a rectangular geometry near its transi- 
tion temperature. The order parameter can behave dif- 
ferently depending on the aspect ratio and size. The 
phase always preserves time-reversal symmetry, except 
for a square with parameters rather far away from 
isotropic weak-coupling theories. Our results give fur- 
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